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Abstract 

We carry out an analytic study of the early-time motion of a quark in a strongly- 
coupled maximally-supersymmetric Yang-Mills plasma, using the AdS/CFT corre- 
spondence. Our approach extracts the first thermal effects as a small perturbation 



D ■ of the known quark dynamics in vacuum, using a double expansion that is valid 



for early times and for (moderately) ultrarelativistic quark velocities. The quark is 
found to lose energy at a rate that differs significantly from the previously derived 
stationary/late-time result: it scales like instead of T^, and is associated with a 
' friction coefficient that is not independent of the quark momentum. Under conditions 

^ ■ representative of the quark-gluon plasma as obtained at RHIC, the early energy loss 

t:J" . rate is a few times smaller than its late-time counterpart. Our analysis additionally 

CN) ! leads to thermally-corrected expressions for the intrinsic energy and momentum of 

the quark, in which the previously discovered limiting velocity of the quark is found 
to appear naturally. 
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1 Introduction and Summary 

1.1 Background: Energy Loss and AdS/CFT 

Heavy ion collision experiments at RHIC and, recently, LHC, have amassed substan- 
tial evidence supporting the detection of the long-sought quark-gluon plasma (QGP), 
a hot and dense phase of deconfined strongly- interacting matter [1]. Energetic partons 
serve as important probes of this thermal medium, and there exists an enormous 'jet- 
quenching' literature dedicated to analyzing the manner in which the plasma damps 
their motion, and is in turn disturbed by their passage [2]. Given the experimental 
indications that the QGP is strongly-coupled, perturbative QCD is believed to be 
inadequate for at least some aspects of the relevant calculations, creating a demand 
for new theoretical tools. In the past five years, interesting steps have been taken 
towards meeting this demand via the AdS/CFT, or, more generally, gauge/gravity 
correspondence [Sill], starting with the seminal works [S1E1[71[S] and continuing with 
a large body of work that has been reviewed in P, [IDl [HI [121 [131 [El- 

The gauge/gravity correspondence is by now well-established as a tool that grants 
access (often analytically) to a large class of strongly-coupled non-Abelian gauge 
theories, via a drastic and surprising rewriting in terms of string-theoretic (frequently 
just supergravity) degrees of freedom living on a curved higher-dimensional geometry. 
The gauge/gravity catalog includes some theories that display characteristic QCD-like 
physics such as confinement and chiral-symmetry breaking, but not, to date, QCD 
itself. Nonetheless, the correspondence has already been very successful in developing 
our intuition on the behavior of strongly-coupled gauge theories, even to the point of 
providing useful suggestions for phenomenological models of the QGP [Ul [ISl [IZl [IB] . 
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In the gauge/gravity context, a strongly-coupled thermal gluon plasma is de- 
scribed by a black hole geometry with special asymptotics, and a quark traversing 
the plasma is dual to a string moving on this background. More precisely, the quark 
corresponds to the tip of a string, whose body codifies the profile of the non-Abelian 
(near and radiation) fields sourced by the quark. As the string endpoint moves, its 
body lags behind it, exerting a drag on the tip that is the gravity-side realization of the 
damping force exerted by the plasma on the quark, as studied in [3 EjQ These works 
focused on the particular case of a plasma in maximally-supersymmetric [Af = 4) 
Yang-Mills theory (MSYM), which is the simplest and best-understood example of 
the correspondence, and is known to furnish a useful (albeit rudimentary) toy model 
of the QGP. MSYM is a conformally- invariant theory (CFT), whose dual gravity de- 
scription at finite temperature involves a black hole that is asymptotically anti-de 
Sitter (AdS). For concreteness, we will hereafter concentrate on this same AdS/CFT 
setup, even though we expect most of the qualitative lessons of our work to be more 
generally 

The analysis of [H [6] led to explicit formulas for the rate of energy and momentum 
radiated by a heavy quark as it moves through an infinite static MSYM plasma. 
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^^XT^^=L= (1) 



dt 2 y/T 



V 



and 



'^■rxr-^J= , (2) 



dt 2 



where A denotes the 't Hooft couplingjfl In the context of weakly-coupled QCD, it is 
known that energy loss at high parton velocity is dominated by the strong-interaction 
analog of bremsstrahlung, i.e., medium- induced gluonic radiation, while for low veloci- 
ties coUisional loss becomes important |2]. In the strongly-coupled MSYM setup made 
available to us by AdS/CFT, the perturbative terminology is no longer adequate. The 
flow of energy/momentum along the body of the trailing string corresponds to en- 
ergy/momentum transported away from the quark via the gluonic field generated 
jointly by the quark and the plasma, in a pattern that has been meticulously studied 
in a large body of work that began with [30], has been reviewed in [3l] and includes 



^It was later shown in [T9l[20] that this traihng string mechanism is in fact much more general, 
and accounts even for the radiative damping expected in vacuum. That is, irrespective of whether 
a spacetime black hole is present or not, the body of the string plays the role of an energy sink, as 
befits its identification as the embodiment of the gluonic degrees of freedom. On the other hand, 
energy loss via the string does turn out to be closely associated with the appearance of a worldsheet 
horizon, as noticed initially in [211 ES] at finite temperature and in [23] (see also [24l [25]) for the 
zero temperature case. 

^Energy loss studies have been carried out via the gauge/gravity correspondence in many other 
theories, including some with properties more akin to those of QCD — see, e.g., [221 [211 references 
therein. Other examples of drag force calculations may be found in the reviews [HI [IHl [HI [HI [131 [E] ■ 

■^Various works have explored the way in which these heavy quark results are modified for heavy 
sources of the gluonic field in color representations other than the fundamental (including the adjoint) 
[28] . as well as for light quarks and gluons [29] . 
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the interesting recent additions [32] . For ease of language, throughout this paper we 
will continue to speak of radiation, even though this concept is not really appropriate 
within an infinite thermal medium. 

We should also stress that, within the framework of O |6], adopted also in this 
paper, the entire energy loss calculation is carried out inside the MSYM theory, 
whose coupling does not run and is taken to be large. For quarks that are extremely 
energetic, the asymptotic freedom of QCD would lead one to expect deviations be- 
tween this scenario and the real-world QGP, and seek a hybrid approach describ- 
ing hard/perturbative emission of a gluon that subsequently propagates through a 
strongly-coupled medium, as advocated in [SI [331 It is debatable, however, 

which of these two scenarios is more appropriate at the not extremely relativistic 
heavy quark energies achieved in RHIC. At the very least, our calculations have a 
direct interpretation in terms of energy loss in the strongly-coupled MSYM plasma, 
which is an interesting theoretical question in its own right. 

1.2 Motivation, outline and main results 

Naturally, expressions ([T]) and ([2]) were derived under a number of simplifying as- 
sumptions. The one that matters most for our purposes is the restriction made by 
the authors of [5], [6] to a configuration where the quark is either moving with con- 
stant (possibly relativistic) velocity as a result of being pulled by an external force 
that precisely balances the drag, or is unforced but moving nonrelativistically and 
about to come to rest. Either of these scenarios requires the interaction between 
the quark and the medium to occur over a considerable period of time. Since the 
actual QGP produced at RHIC or LHC has a finite temporal and spatial extent, it 
is questionable whether a heavy quark propagating through it will have enough time 
to reach a quasi- stationary configuration, or to be nearly stopped (in which case it 
would anyhow not lead to an experimentally distinguishable jet). 

The actual energy/momentum loss might thus be expected to differ from ([T]) 
and ([2]) in a situation where the quark moving through the plasma is accelerating, 
or in the initial period following its production within the thermal medium. This 
point has been emphasized from the phenomenological perspective in [381 EH]. The 
estimates there are based on perturbative calculations, so it is interesting to inquire 
into this issue in the strongly-coupled systems available to us through the AdS/CFT 
corresp ondence . 

Now, the restriction in [5l [6] to the stationary or asymptotic cases was of course 
implemented to gain analytic control on the problem of energy loss. The authors 

^The approach of O |33] involves the extraction from a hghthkc Wilson loop of a transport 
('jet-quenching') parameter which characterizes the medium and whose value at strong coupling 
curiously differs from the customary definition of q as the average transverse momentum picked up 
per unit distance traveled [HI [52] . Some doubts about the procedure proposed in [51 [33] to calculate 
q on the gravity side of the correspondence had been raised by, e.g., |34l [35l [36], but these concerns 
appear to have been resolved by a recent correction to the procedure [37], which happens to leave 
the result unchanged. 
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of [IQ] were able to slightly extend the stationary result to the case of a slowly de- 
celerating quark, working in an expansion in powers of \/XT/m. Away from these 
regimes, it is difficult to follow the evolution of the quark in the thermal plasma, or 
equivalently, of the string on the AdS black hole geometry. One possibility is to re- 
sort to numerical analysis, a strategy adopted in [23]. Among other things, that work 
studied the case where a quark that is initially static within the plasma is accelerated 
by an external force over a finite period of time and is thereafter released. It was 
found that under such conditions, and for values of the mass in the neighborhood 
of the charm quark, there exists a period after release where the quark dissipates 
energy at a rate that is substantially smaller than the stationary/asymptotic result 
([1]) obtained in [5l |6]. In addition, the rate of energy loss in the initial stage where the 
quark is externally forced was found to be almost completely accounted for by a gen- 
eralized Lienard formula describing radiation in vacuum (derived in that same paper 
and in [191 ISO])- Unfortunately, the numerical integration degraded rather quickly, 
as a result of which the investigation of [23] was limited to intervals that are more 
than an order of magnitude below the experimental timescale of the plasma (typically 
^breakdown ~ O.O/ttT ~ 0.3 fm/c), and to rather small quark velocities at the time of 
release (for the most part, ^release < 0.1). 

Both to try to overcome these numerical limitations, and for the additional intu- 
ition it would provide, it is clearly desirable to gain some measure of analytic control 
over the problem of early-time quark damping in a strongly-coupled plasma, which is 
what we set out to accomplish in the present paper. Our approach will be to study 
thermal effects as a small perturbation on the zero-temperature evolution, which is 
reasonable on general physical grounds during the initial stage of motion through the 
plasma, and is moreover supported by the numerical results of [23]. For simplicity, we 
will restrict attention to the case where the quark moves solely along one direction, 
denoted by x, as this is enough to examine the damping we are after. 

After setting up our language and notation in Section |2] (which can be safely 
skipped by the cognoscenti), we begin to develop our story by establishing the zero 
temperature framework over which we will perturb. In Section [STTl we recall the string 
embedding ( ITT]) , argued to be an extremum of the Nambu-Goto action in [H], and 
explicitly verify that it indeed solves the corresponding nonlinear equation of motion. 
This provides the string configuration relevant for an arbitrary timelike motion of 
the heavy quark in vacuum, under the condition that waves on the string be purely 
retarded, to capture the gluonic fields causally set up by the quark. In Section |3^ we 
then present a streamlined version of the procedure of [231 [IHl 120] that uses the string 
profile ( irTj) to turn the standard string boundary condition f|T^ into a dynamical 
equation for the quark in vacuum. The result, equation flHT]) . is a generalized version of 
the classic Lorentz-Dirac equation [THl [20] that is physically sensible, and incorporates 
expressions for the intrinsic momentum of the quark (!32|) and the rate (133|) at which 
it radiates momentum. Equivalently, equation (fT6|) incorporates expressions for the 
intrinsic quark energy (!36|) and rate of energy radiation (1371) . As found in [231 UHl |20] , 
for a quark with finite mass and therefore finite size, these expressions differ from the 
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naive ones by corrections that depend on the external force and the quark Compton 
wavelength. In the quark equation of motion we can solve for the acceleration to 
arrive at fl40l) . which shows that, as expected (and unlike in the standard Lorentz- 
Dirac equation), an unforced quark in the vacuum of MSYM must move at constant 
velocity. This is of course the property that we expect to be modified upon the 
introduction of a thermal medium. 

In Section 14.11 the finite-temperature string equation of motion fH^ is considered 
for early times in the evolution, when the disturbed portion of the string remains 
far from the black hole horizon, as indicated in (H3|) . and with the embedding (!52|) 
deviating only slightly from its zero-temperature counterpart (fT7|) . to arrive at the 
linearized equation (146|) . This still contains a very complicated dependence on the 
quark trajectory, so, motivated by the experimental scenario, we further restrict at- 
tention to quarks traveling rapidly through the plasma. In short, we work in the 
early-time and moderately ultrarelativistic regime f H9|) . which essentially amounts to 



TT^T^ci^ < 1 - < 1 , (3) 

where d ~ t/7 is the size of the growing region of the gluonic field that has been 
disturbed by the motion of the quark. Under these conditions we are ultimately led 
to the thermally-corrected embedding fl6T]) . 

We then go on in Section U2]to deduce the way in which the first thermal correction 
to the embedding affects the evolution of the quark, focusing attention on the case 
of prime phenomenological relevance where the quark is not externally forced. By 
imitating the procedure used in the zero-temperature setting, we derive an equation 
of motion for the quark in the thermal medium. Along the way, we find that our 
framework naturally makes contact in (|74l) with the previously known limiting velocity 
Urn < 1 (defined in ( 17T]) ) for a quark immersed in the strongly-coupled plasma [121 
2Il|221i2llH[3ni[2SlllSlllS]. The thermal equation of motion that results from 
our analysis can be written in either of the alternative forms flHT]) . f l38|) or flMjl . with 
the identification of f l85|) and f l83|l as the thermally corrected intrinsic energy and 
momentum of the quark (which are seen to diverge as v — v^)-, and of 



dt 16 m? 1 — {v/vmy 

and 

dR., TT 1 



(4) 



(5) 



dt 16 1 — {v /VraY 

as the rates at which the quark loses energy and momentum. These are then our 
main results. 

Section |5] discusses the possible phenomenological implications of our findings. 
After verifying that our approximations are well-suited to the typical energetic heavy 
quarks obtained at RHIC, we contrast our early-time results (|4j)-(l5]) against the late- 
time expressions ([I])-(l2]) obtained in [HIE]. In ( 19T|) we see that, for parameter values 
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appropriate to RHIC, the former can be up to six times smaller than the latter. This 
comparison is consistent with the (more hmited) numerical results of [23]. The func- 
tional form of our results is also of interest. From our derivation it becomes clear that 
the T^-dependence seen in dl])-© (and anticipated in an estimate based on satura- 
tion physics [21]) follows naturally from the first thermal corrections in the metric, 
and via dimensional analysis, this fixes the accompanying dependence on the quark 
mass, which in turn explains the peculiar power of the 't Hooft coupling. Another 
prominent feature of our formulas is their velocity-dependence: whereas the late-time 
friction coefficient ( 193|) deduced from (II])- ([2]) is momentum- independent, its early-time 
counterpart (l94l) depends strongly on momentum. This feature should constitute an 
interesting experimental signature when used as input for phenomenological models, 

as in [TsmniiniiiH]. 

As envisioned from the outset, our investigation is limited to the initial stage where 
the quark is only starting to feel the effects of the plasma. In fl95|) we learn that our 
results should be reliable for times that are not too close to the lifetime of the real- 
world QGP. The temporal range of validity of our work is an improvement over the 
one accessible in the numerical work of [23], but still suggests that the experimental 
results should effectively arise from some sort of average or interpolation between the 
early and late energy loss rates (j4]) and ([1]). In this connection it should also be kept 
in mind that a separate numerical study in [23] (building on [5]) of a quark that is 
created within the plasma back-to-back with its corresponding antiquark found the 
late-time energy loss rates of [HI [6] to be relevant essentially as soon as the quark 
and antiquark separated beyond the velocity-dependent screening length determined 
independently in [171 El] (see also [IH]) and [23] itself. The way in which these findings 
interface with ours is also discussed in Section |5l The upshot is that, in the more 
realistic situation where the presence of the accompanying antiquark is taken into 
consideration, our results would be most relevant in the stage before the separating 
quark and antiquark lose contact with one another. This is natural, because screening 
emerges only when thermal effects are substantial. For a quark ploughing through a 
significant portion of the plasma, then, the late-time formulas ([I|)-([2|) actually stand 
a chance of controlling the larger portion of the evolution. 

Although we have succeeded in gaining some analytic control on parton energy 
loss in a plasma within a novel temporal regime, clearly more work will be needed to 
shed additional light on the overall problem and to remove the various idealizations 
adopted in the existing calculations. Based on past experience, it seems likely that 
the gauge/gravity correspondence will continue to be a useful guide into the strong- 
coupling aspects of this problem. 

2 Preliminaries: Quarks as Strings 

Maximally supersymmetric (i.e.. A/" = 4) Yang- Mills (MSYM) is a conformally invari- 
ant theory comprising a gauge field, 6 real scalar fields and 4 Weyl fermions, all in 
the adjoint representation of the gauge group. According to the AdS/CFT correspon- 
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dence [3], large- A^^^ strongly-coupled SU{Nc) MSYM on (3+l)-dimensional Minkowski 
spacetime, with 't Hooft coupling A = Qym^c and temperature T, is equivalent to 
Type IIB string theory living on the (planar Schwarzschild-AdS)5 x geometry 

ds^ = Gmndx'^dx'' = — -hdt^ + dx^ + —]+ R^d^l , (6) 

\ h J 

(with a constant dilaton and N^. units of Ramond-Ramond five-form fiux through the 
five-sphere,) where Is denotes the string length. The radial direction z is mapped 
holographically into a variable length scale in the gauge theory |19] . The coordinates 

= {t, x) are parallel to the AdS boundary z = and are directly identified with the 
gauge theory spacetime coordinates, while the five-sphere coordinates are associated 
with the global SU (4) internal (R-) symmetry of MSYM. 

The state of IIB string theory described by the pure AdSs x geometry given 
by the z^ ^ oo limit of (|6]) corresponds to the vacuum of the MSYM theory, and the 
closed string sector describing (small or large) fiuctuations on top of it fully captures 
the gluonic (+ adjoint scalar and fermionic) physics. To this closed string sector 
belongs in particular the Zh < oo geometry ([6]), which describes a black hole (strictly, 
black brane) with an event horizon at 2; = 2;^ and Hawking temperature T, and is 
dual to a MSYM thermal ensemble at the same temperature. This then is the string 
theory implementation of our desired plasma of gluons (-1- scalars and fermions). 

To this thermal MSYM setup we can add Nf hypermultiplets (each composed 
of a Dirac fermion and 2 complex scalars) in the fundamental representation of the 
SU{Nc) gauge group, breaking the supersymmetry down to A/" = 2. These are the 
degrees of freedom that we will refer to as 'quarks,' even though they include both 
spin 1/2 and spin fields. In string theory language, adding these hypermultiplets 
corresponds [50] to the introduction of an open string sector associated with a stack of 
Nf D7-branes in the geometry ([6]). For Nf <^ Nc, the backreaction of the D7-branes 
on the geometry can be neglected; in the gauge theory this corresponds to working 
in a 'quenched' approximation that ignores quark loops (as well as the positive beta 
function they would generate). 

The D7-branes cover the four gauge theory directions t, x, and extend along the 
radial AdS direction up from the boundary at z = to a position where they 'end' 
(meaning that the S'^ C that they are wrapped on shrinks down to zero size), 
whose location z = Zm is related to the mass m of the quark through [5] 



1 27rm 



V'A 




128 \ 2m I 2m 



(7) 



We are interested in quarks that are heavy in the sense that m > y/\T/2, so we will 
restrict attention to D7-brane embeddings with 2^ -C 2^, and in so doing will only 
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need the leading order approximation to (j7]), 

Zm — „ 

zvrm 



(8) 



An isolated quark is dual to an open string that extends radially from the D7- 
branes to the black hole horizon aX z = z^. We take the string to lie, consistently 
with its equations of motion, at the 'North Pole' on the (the point where the 
S'^ C that the D7-branes are wrapped on collapses to zero size), so the angular 
components of the metric will not play any role in our work. This is natural for 
eventual application to the real-world QGP, because QCD has no symmetry on a par 
with the SU{4) internal symmetry of MSYM that is dual to the five-sphere isometry 
group. The lower endpoint of our string will then necessarily be located aX z = Zm- 

As is customary in string theory, we will describe the dynamics of our string in 
first-quantized language, and since we take it to be heavy, we are allowed to treat 
it semiclassically. In gauge theory language, then, we are coupling a first-quantized 
quark to the gluonic (+ other MSYM) field(s), and carrying out the full thermal 
path integral over the strongly-coupled field(s) (the result of which is codified by 
the Schwarzschild-AdS spacetime), but treating the path integral over the quark 
trajectory x{t) in a saddle-point approximation^ 

It is useful to keep in mind that in the nonperturbative context provided to us 
by the AdS/CFT correspondence, a quark with finite mass > 0) is automatically 
not 'bare' but 'composite' or 'dressed'. This can be inferred, for instance, from the 
expectation value of the gluonic field surrounding a static quark [57], or from the 
deformed nature of the quark's dispersion relation [23l [T9l [20] . As seen in those 
calculations, the characteristic thickness of the 'gluonic cloud' surrounding the quark 
is given by z^, which is thus understood to be the analog of the Compton wavelength 
for our non-Abelian source. 

The string dynamics follows as usual from the Nambu-Goto action 

'S'ng = J a/- det gab = J c^V^ng , (9) 

where gab = daX"^dbX^Gmn{X) {ci,b = 0,1) denotes the induced metric on the 
worldsheet. For convenience we choose to work in the static gauge ct° = t, = z. 

By the AdS/CFT dictionary, it is the endpoint of the string on the D7-branes 
that directly corresponds to the quark, and so the latter's trajectory can be read off 
as 

x{t)=X{t,z^) . (10) 

In our analysis below we will have need to refer to the velocity, acceleration and jerk 
of the quark, which will be denoted a.s v = dx/dt, a = dv/dt, j = da/dt. Since the 



^Starting with (HH [52], various works have explored quantum fluctuations of the quark about 
certain types of average trajectories — see, e.g., [53l[54l[55] and references therein. These fluctuations 
are suppressed by a factor of I/v'A, and would have to be complemented with contributions from 
the path integral of the type explored in a somewhat different context in [56] . 
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damping force we are after will be exerted by the plasma against the quark's direction 
of motion, for simplicity we will consider motion and deformation of the string purely 
along a single direction x = x^. Extension of our entire study to arbitrary three- 
dimensional motions does not require any new conceptual ingredients (and at zero 
temperature has already been carried out in [T9l [20] ). but is algebraically much more 
complicated (especially since at finite-temperature one can no longer rely on Lorentz 
invariance as in [121 |2U] ) . 

Under the simplifying conditions described in the two previous paragraphs, the 
Nambu-Goto Lagrangian boils down to 




£ng = + - ^ , (111 
and the non-zero canonical momentum densities IIJJ = dCj<sG/d{daX'^) are given by 



hX'^ + 1 



z^sJl + hX'^ 


X2 

h 


X 




z^h^Jl + hX'^ - 


h 


hXX' 




z^^l + hX'^ - 


) 

h 


hX' 





III = , , (12) 



z^^ll + hX''-^ 



where of course X = dtX{t, z), X' = dzX{t, z). Notice that, due to our normalization 
of £ng, the n° must be multiplied by i?^/27r/^ = -\/A/27r to obtain the physical energy 
and momentum densities. 

We can exert an external force F on the quark/string endpoint by turning on an 
electric field = Fi on the D7-branes, which contributes to the string action the 
term 

Sf = J dtA^{X{t,Zm))dtX^{t,zJ . 

Variation of the string action Sng + 'S'f implies the Nambu-Goto equations of motion 
(or, equivalently, current conservation conditions) 

d,ul + a,n^ = (13) 

for all interior points of the string, plus the standard boundary condition 

n^WI.=.™ = ^-^.W vt, (14) 

where J-"^ = —F^^dtx'^ = {—F-v, F) is the Lorentz four-force. In the following sections 
we will solve for the dynamics of this system to extract the thermal damping rate. 
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3 Zero Temperature 



In this section we will examine the quark dynamics in the MSYM vacuum, to lay the 
groundwork for the finite-temperature case that is our real interest. In so doing, we 
will rederive a subset of the results of [HI [231 IIHl 120] in a form that is better suited 
to our purposes. 



3.1 String embedding 

At zero temperature {z^ — )■ oo) the background reduces to pure AdS, and the 
string equation motion f lT3|) reads 

\ a / VI \ 

= , (15) 





which can be reprocessed to 

X-X" + XX'^ + X"X^ - 2XX'X' + + - = . (16) 

z z z 

To understand the behavior of the quark in the MSYM vacuum, our task is then 
to solve equation f fT6|) subject to the forcing boundary condition flT^ . and read off 
the quark trajectory from the result using ( ITOi) . Equivalently, we can take a specific 
quark trajectory and interpret it through ( ITOi) as a boundary condition for the string, 
solve for the corresponding string embedding, and in the end read off from ( IT4l) the 
forced required to guide the quark along the chosen trajectory. For the time being 
we will adopt this second perspective and, moreover, will choose to parametrize our 
desired string embedding in terms of the data for a quark that is infinitely massive, 
thus setting Zm = until further notice. In |?T], Mikhailov used the symmetries of 
the AdS geometry to argue that, for an arbitrary time-like quark trajectory x(t), the 
embedding 

X{U,z) = ^itr) + ^fi^^^, (17) 

t(^trp^ tiy I 



which manifestly meets the condition (fTOj) . also extremizes the Nambu-Goto action, 
and must therefore be the solution we are seeking. 

From the structure of (1171) we see that the behavior of the string at a given time 
t and radial depth z is completely determined by the behavior of the quark/string 
endpoint at the earlier, retarded time tr- The definition of tr implicit in f lTTj) can be 
shown to follow from projecting back to the AdS boundary along a curve that is null 
on the string worldsheet [H], in analogy with the Lienard-Wiechert story in classical 
electrodynamics. The embedding ffTT]) thus describes a wave on the string that is 
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purely outgoing: it is generated at the string endpoint and then ascends along the 
body of the string, moving into the AdS bulk. This is why, among the infinite number 
of extremal string embeddings that satisfy fllOl) . the profile fll7l) is the one that is of 
immediate physical interest for us: it is dual to a MSYM configuration where waves 
in the gluonic field move out from the quark to infinity. 

We will now show explicitly that f lT7|) indeed solves the equation of motion f|T6|) . 
Noting that 



dt 
dX 



dtr 
dtr 



vaz 



az 



+ 1 



+ 



dz 



v'^az 



^2)3/2 



vdz 



where from now on it is understood that the velocity, acceleration and jerk of the 
quark (in the notation introduced after (fTOj) ) are evaluated at tr, we can deduce after 
some algebra that 



X 



fdX_ 



'dt 



fdX_ 




dx_ 

dz 



az + v{l - v'^f/'^ 

vaz + (1 — t>2)3/2 ' 

az{l - t;2)i/2 

vaz 



(19) 
(20) 



1 - t;2)3/2 ' 

-a3^2(l - t;2)3/2 + 3m2z(l - t;2)3 + jz{l - v^f + a(l - v^fl^ 

[uaz + (1 - t;2)3/2]3 

-a3^2(l _ .^;2)l/2 ^ ^^2^^i _ ^2)2 ^ ^^^-^ _ ^,2)3 _ _ ^2)7/2 



a^zHl 



[vaz + (1 - i;2)3/2]3 

;2) - 2va''z{l - v^f'^ - jz{l - v^l^ 



\vaz 



2)3/2j; 



(21) 



Substituting ( |T9l) . ( l20l) and ( 1211) into the equation of motion (|T6l) . it is easy to verify 
that it is in fact satisfied, thus providing an independent check of the extremization 
argument of |41j . 



3.2 Radiation damping 

When the quark/endpoint moves at constant velocity v, (IT7|1 describes an upright 
(i.e., purely radial) string translating uniformly, as one could have anticipated by 
Lorentz invariance. For any accelerated motion, an event horizon develops on the 
string worldsheet [23] , associated with the expected fact that the quark emits gluonic 
radiation and is thus subjected to a corresponding damping force. In [T9l [20] it was 
shown, for general three-dimensional motions, that this physics is correctly captured 
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by the forced boundary condition (fT^ . Indeed, through use of Mikhailov's embedding, 
this condition can be rewritten purely in terms of the time evolution of the endpoint, 
and so can be interpreted as an equation of motion for the quark, incorporating 
Lorentz-covariant formulas for its rate of radiation and intrinsic four-momentum. We 
will here reobtain this information via a streamlined procedure adapted to our one- 
dimensional setting. 

For the radiation damping to be noticeable at the level of the quark equation 
of motion, we need the quark to have a finite mass, so from now on we allow again 
Zm> ^ (which, as we emphasized in Section 2, inevitably means that our non-Abelian 
source is no longer pointlike but has size Zm)- The novelty is then that the boundary 
condition ( ITO|) is now imposed not at the AdS boundary but at the radial position 
z = Zm where the string endpoint lies. As before, this condition by itself does not 
pick out a unique string embedding. Just like we discussed in the previous subsection 
for the infinitely massive case, we additionally require the solution to be 'retarded' or 
'purely outgoing', in order to focus on the gluonic field causally set up by the quark. 
As in [23| [T9| 120] . we can inherit this structure by truncating a suitably selected 
retarded Mikhailov solution. The solutions of interest to us can thus be regarded as 
the z > Zm portions of the embeddings fll7p . which are parametrized by data at the 
AdS boundary 2; = 0. From this point on we will use tildes to label these (now merely 
auxiliary) data, and distinguish them from the actual physical quantities (velocity, 
time, etc.) associated with the endpoint /quark aX z = Zm, which will be denoted 
without tildes. 

In this notation, ( IT71) reads 

X{ir,z) = x+ , (22) 

V 1 — w 

t(^tf^^ — t^" I 



and expressions 0191) . fl20p and fl2ip continue to hold if their right-hand sides are 
similarly decorated with tildes. Using this information in (IT^ we can deduce that 

^ z{l - {;2)3/2 ' ^'^'^^ 
which when evaluated ai z = Zm implies, with the aid of ([H]), that 

—F = f . (24) 



zm{i-v^r/' VX 

By evaluating (fT9|) at the endpoint we can also infer that 

azm + v{l-v'^)'^/^ 

V = — 

vaZm + (1 — f 2)'^/2 



(25) 
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which we can invert in combination with (1241) to learn that 

^2)3/2(1 4^ ^2)3/2 



(26) 



Having aheady incorporated the dynamical information contained in (fT4|) and 
(fT7|) . the equation controlling the evolution of the quark now follows simply by sub- 
stituting ( 126|) into the simple kinematical statement 



dtr 



Reading off from (ITTl) that 



tj. — t(tf., Zyu) -\- , (2J 

V 1 — 



we can deduce with the aid of ( l26l) that 

dtr _ (1 - v^/^ + ^az^ _ 1 - z^f ^ 

dir (1-^2)3/2 l-Z^t^f 

Together with (12^ . this turns ( 127|) into 



(29) 



or, equivalently, 



(31) 



which is the quark equation of motion we were after. As explained in [13 |2U], it 

has a transparent and pleasant interpretation in gauge theory language: recognizing 

rri'yv — -^^^F 

Pg = ^""^ = (32) 



A 



as the intrinsic momentum of the quark and 



_ P2 

4 -f^ 



dt 27r m2 I 1 - 



(33) 
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as the rate at which momentum is carried away from it by the emitted gluonic radi- 
ation (and consequently, as the negative of the radiative damping force), (131 p is seen 
to amount merely to the statement of momentum conservation. 



dP _ dpg rfPrad 



(34) 



dt dt dt 

The total momentum supplied to the system by the external force must either increase 
the intrinsic momentum of the quark or be radiated away into the gluonic field. The 
fact that this eminently physical gauge-theoretic statement has arisen so straightfor- 
wardly from the standard string dynamics is a nice illustration of the power of the 
AdS / CFT correspondence. 

Needless to say, the total energy is also conserved, and the corresponding state- 
ment can be obtained by applying the above argument to d^/dir = ^^vd instead of 
(1271). yielding 



m7 




Identifying 



m7 — ^^'yvF , ^ 

Eg = ^^"^ ' - (36) 



as the intrinsic energy of the quark and 

dE,^ _ VXF' f 
dt 2vrmMl-^F2 



(37) 



as the rate at which it radiates energy, (l35ll embodies the energy conservation state- 
ment 

dE dEg dE^ad T? /oo\ 

—— = —-^ H ; — = vF . (38) 

dt dt dt ^ ^ 

The Lorentz-covariant generalization of equations fl3T|) and fl35l) . for arbitrary 
three-dimensional motions of the quark, can be found in [T^ |2U]. As emphasized 
there, the result can be understood to be an extension of the well-known Lorentz- 
Dirac equation of classical electrodynamics [58] to the strongly-coupled non-Abelian 
setting. Unlike its Lorentz-Dirac counterpart, the AdS/CFT extension is found to be 
non-linear and to have no pathological solutions. 

We note in passing that (1261) and (l28l) allow us to eliminate all auxiliary quantities 
and rewrite the string embedding (ITTl) purely in terms of physical quark/endpoint 
data, 

t{t.,z) = t.+^-^-^^ntzM. 
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The non-trivial form of the quark dispersion relation embodied in (132|) and (!36l) 
is a direct consequence of its extended, non-pointlike, nature, which we took care 
to emphasize in Section [21 Indeed, the characteristic length scale appearing in the 
quark dressing seen in fl32|) is none other than Zm, which we know to play the role 
of Compton wavelength. The AdS / CFT correspondence is thus teaching us precisely 
how the gluonic cloud surrounding the quark is deformed upon the application of an 
external force. 

For F = 0, f l5^ reduces to the usual dispersion relation pg = m'-fv (and similarly 
Eg = Twy). The same is true in the infinite-mass/pointlike limit Zm ^ 0. A conse- 
quence of this last statement that we would like to stress as an important lesson of the 
above procedure is that, in terms of the auxiliary/tilde variables, the intrinsic quark 
momentum and energy always take the usual form, pg = rri'yv (as seen in ( 1271) ) and 
Eg = m7, and the non-trivial dressing present in 0321) and (1361) arises only upon using 
fl26l) to rewrite this in terms of the physical/nontilde variables. In the next section 
we will learn that this continues to hold true even in the finite-temperature setting. 

For later use it is also useful to recall that in [HI [23] , a time- integrated version 
of the momentum/energy split seen in (I34p and (1351) was shown to be achieved on 
the string theory side directly at the level of the total momentum/energy carried by 
the string worldsheet. In more detail, [H] obtained the infinite-mass limit of (l33l) 
and (1371) (which, surprisingly, has the same functional form as the Lienard formula 
of classical electrodynamics) by using (1T71) in the integral of the momentum/energy 
densities 11^ and 11* given in (lT2l) . These calculations disregarded a total derivative 
term, which was later shown in [23] to correspond to the infinite- mass limit of (l32l) . 
( l36l) . The same reference extended the analysis to the case of finite mass. Of this 
story, the element that we most want to emphasize is that, at the worldsheet level, the 
intrinsic quark momentum and energy (l32l) and (l36l) are obtained as surface terms, 
which are naturally evaluated at the string endpoint. 

Our one-dimensional equation of motion ( 1311) or (1351) can be rewritten in the form 



from which it is clear that the quark will accelerate (and thereby radiate) only if 
externally forced. This is precisely as one would expect for a particle in vacuum, but 
of course the situation will be different when we introduce a thermal medium in the 
next section. 

4 Finite Temperature 

We will now extend the results of the previous section to the case where the quark 
moves inside a thermal MSYM plasma. In Section 14.11 we will obtain the first finite- 
temperature correction to the string embedding, and then in Section 14.21 we will use 
this information to infer the rate at which the quark is damped by the plasma. 



z^ra-v^/' 



+ 



(40) 



a = 
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4.1 First thermal correction to string embedding 

When we turn the temperature back on we have Zh = l/vrT < oo, and the momenta 
f lT2|) imply that the string equation of motion f fT3|) now takes the form 



X \ d hX 



, (41) 



z'^h^l + hX''^ - X'^/h) \z'^\ll + hX''^ - X^/h^ 

or, equivalently, 



= X - h^X" + hXX'"^ + hX"X'^ - 2hXX'X' + 
_ hh'X' - ^h^h'X''' + Ih'X'X^ . 



2h^X' 2h^X'^ 2hX'X^ 
~^^^z i (42) 



It would be highly desirable to obtain from this finite-temperature equation a 
family of exact solutions analogous to f|T7|l . but unfortunately, we have not been 
able to accomplish this. Our aim will be instead to derive the first corrections to the 
string/quark dynamics induced by the presence of the thermal medium. As explained 
in the Introduction, when a heavy quark (m ^ \/AT/2) is created within the MSYM 
plasma we expect its initial evolution to be quite close to that in vacuum, with the 
thermal effects only starting to be felt gradually, as the gluonic field set up by the 
quark begins to interact with the medium and thereby induce a damping effect on 
the quark. In the string theory description, this corresponds to the fact that (for 
Zm/zh ^ 1) a wave propagating upward from the string endpoint is initially far from 
the Schwarzschild-AdS horizon, and so evolves essentially as in pure AdS, with small 
corrections arising gradually as the wave moves deeper into the bulk. As long as we 
restrict attention to the early-time evolution where the disturbed portion of the string 
remains far from the black hole horizon, 

TT^T^z^ < 1 , (43) 

the form of the metric (jH]) enables us to regard the string profile as a small pertur- 
bation about the zero-temperature solution f[T7|) . By the UV/IR connection [IH], z 
corresponds to a length scale d = z on the MSYM side, so the condition fH5]) states 
that the region of the gluonic field that is disturbed by the motion of the quark re- 
mains small in units of the inverse temperature of the plasma. With this restriction 
in mind, we thus expand the exact equation of motion ( lT3l) in powers of the small 
parameter ( H3l) and keep only the leading non-trivial correction: 



= \X-X" + XX" + X"X^ - 2XX'X' + ^ + ^- 

z z z 

+ TT^T^z' I 2X" - - X"X^ + 2XX'X' - - ^ 1 + Oin'T^'z' 



(44) 
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In parallel with this, we write the string embedding as a perturbation of Mikhailov's 
vacuum profile (1171) . which we hereby denote with a "0" subindex, 

X{t,z)=Xo{t,z)+Y{t,z) , (45) 

and where it is understood that |y/Xo| <^ 1 throughout the region of interest. Plug- 
ging this ansatz into (jSj), and keeping only terms up to first order in the perturbation 
Y, we obtain the linearized equation of motion 

AY + BY' + CY + DY" + EY' + F = Q (46) 



where 

A = 2XqXo — 2XqXq 



B — 2XqXq — 2XqXq 

c = x',' + i, 

'2 



Z 

2 ^X'^ 2X1 



D = X^-1, (47) 
E = -2XoX'q , 

F = tt^tV ( 2X^' - XoX^' - X^'X^ + 2XoX^X^ - - 

\ z z 

The coefficients obtained after substituting (fT9l) . (120|) and (!2T!) into ( H7|) have a 
very complicated dependence on v, a, j and z, so it is challenging to find a full 
solution of even the linearized equation (l46ll . We therefore need to make some sort 
of additional simplification. Guided again by phenomenological considerations, we 
choose to restrict attention to a quark that is ultra-relativistic. More specifically, we 
assume that 

a = Vl-v^ < 1 , (48) 

which we will be able below to reexpress directly in terms of the physical quark 
velocity v. We can then proceed to expand ( 1471) in a power series in this additional 
small parameter, using v = \/l — ~ 1 — ^«^, v"^ = 1 — o? ^ -y^ ~ 1 — etc. As 
we will find shortly, the expansion of the embedding will actually end up proceeding 
in powers of o? . Notice that, to the extent that we wish to consider successive 
corrections in o? while continuing to disregard the higher-order thermal corrections, 
we are restricted to the regime where 

tt^TS^ < 1 - < 1 , (49) 

which as we will see in the next section is in fact consistent with the phenomenological 
setup at RHIC. 

Now, in performing the ultrarelativistic expansion of f l47|) . given the form of the 
denominators in f|T9|) . fl20|) . fl2T|) . it becomes necessary to establish whether the quan- 
tity az, which we expect to be small for a thermally-induced acceleration, is smaller 
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than, of the same order as, or larger than the small quantity (1 — f ^)^/^. We have 
explored all three possibilities and found that the only one that is self-consistent is 

« 1 . (50) 



(1 - 52)3/2 

Notice that this constraint is additionally consistent with the facts that (i) we have 
already placed an upper bound on z through (ii) we need the acceleration not 
to be to large in order to remain within the ultrarelativistic regime ( HHj) : (iii) for 
reasons to be explained in Section [5], we will need the disturbance on the string to 
stay below the finite-temperature generalization of the dynamical worldsheet horizon 
that at zero-temperature lies close to a stationary limit curve [23] located at Zcrgo — 
(1 — v^f'I'^la. An additional observation similar in spirit to this last remark is that, in 
the constant-velocity thermal embedding of [5l[6], = {\—v^^^l^Zh marks the location 
of a stationary limit curve that is simultaneously a worldsheet horizon [211 [22], ^ind 
our condition (H9|l restricts us to radial locations below this. 

Combining these assumptions we find that the coefficients (1471) reduce to 

A = 1 , 
5 = a% 

C = (51) 

z 

D = — 

z 

E = 2a 



It is natural then to propose a series expansion for the perturbation on the string 
embedding, 

Y{tr,z) = J2Mzh' , (52) 

k 

where the functions fk{z) could generically also depend on tr through a, j, etc. Upon 
taking the corresponding derivatives in ( 14^ and matching coefficients of a given order 
in a, it is easy to see that the terms with k < —1 lead to homogeneous differential 
equations for f^, meaning that these terms do not contain any information about the 
thermal medium, and can thus be set to zero. (They must in fact be set to zero if 
we are to preserve the purely-outgoing structure of the string profile.) The leading 
order term in fl52l) is then of order 1/a, which as seen in (1221) is the same velocity 
dependence present in the unperturbed solution Xq. 

The k = —1 coefficient is found to satisfy the differential equation 



whose general solution is 



2/:i(^)+7rW = , (53) 



tt'^T'^z^ 

/-i(^) = 1^ + ^- (54) 
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Demanding that Y{t,z = 0) = 0, so that x(t) continues to represent the trajectory 
of the (auxihary or physical, depending on whether or not = 0) string endpoint, 
we are forced to set c = 0, and conclude that, to this order in our small-temperature 
and large-velocity expansion, 

Y{t,^z) = - . (55) 

7 .Z 



lOv/T 



V 



Notice that |y/Xo| is consistently small, of order f l43|) . as expected from the size of 
the inhomogeneous term in the equation P6|) . 

Starting from fl55|) . we can recursively obtain higher-order corrections to the per- 
turbation (!52|) . The /c = coefficient is found to obey 

/oW = , (56) 



implying that 
At level /c = 1 we find 

whose solution is 



/o(^) = . (57) 
2/{(z) + ^7rW = 0, (58) 

A(^) - . (59) 

The next coefficient, f^-, is again found to vanish, so, up to terms of order o? ^ we learn 
that 

^(*-^>^-WT=1J S • 

and therefore 

, vz -n^T^z^ Zi^^T^z^^fY^ , ^ 

X(tr,z)=x^ , , . 61 

One can continue this expansion to higher order in a, and upon doing so one finds 
that terms involving a and j begin to appear in the equation of motion, indicating 
the need for a more general ansatz. We will not pursue this further here, however, as 
our main interest lies in determining the leading thermal damping, which as we will 
now show is available to us already from the first T-dependent term in (16T|) . 

4.2 Early-time thermal damping 

We now use the thermally-corrected string embedding derived in the previous subsec- 
tion to learn about the way in which the MSYM plasma begins to damp the quark. 
Including only the leading nontrivial correction in the small quantities fHSjl and fH8|) . 
we have learned that 

X ~ vz 'k'^T'^z^ 
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From this it follows that 



Mu, .) . (f ) ^ ^"-«'>!'t;-(^;|r-''^) (63) 

\ ot J vaz + {l — v^y'^ 

and 

^ ' \dz)^ y~az^l^ + {I - v'^Y 

which evaluated at the string endpoint yield expressions for the physical quark/endpoint 
velocity v = X(tr, Zm) and the endpoint slope s = X'{tr, Zm)- 

Of particular phenomenological relevance is the case where the quark is not exter- 
nally forced, where we can see from f[T^ and f lT^ that the string satisfies the familiar 
free-end (Neumann) boundary condition s = 0. We can then invert (!63|) and (!64l) 
evaluated aX z = Zm to find that 



V + |vr^T^4 (65) 



and 



10(1 - v^) - Qn^T^zf^v - Itt^T^z^ ' ^^^^ 

It is now useful to remember that, for arbitrary z^, our calculational framework 
interprets the physical string as merely the [z > Zm) portion of an auxiliary string that 
reaches all the way down to the Schwarzschild-AdS boundary. In this interpretation, 
V denotes the velocity of the (fictitious) endpoint at 2; = 0, while v represents the 
velocity of that same embedding at radial depth z = Zm- This implies that we can 
move along the physical string by fixing a value of v (to select a particular embedding) 
and considering Zm as the free parameter z that is greater than or equal to the value 
([7]) determined by the quark mass. We see in fl65|) that, in this reading, the velocity v 
decreases increased, which tells us that, due to the presence of the black hole, 

the horizontal velocity of the string is smaller for points deeper into the Schwarzschild- 
AdS bulk. In other words, whereas at zero temperature an externally unforced string 
would be upright, we find here that the thermal effects make the string lean back, 
and therefore exert a damping force on its endpoint, just as one would expect, and 
as is visible (for the auxiliary endpoint) also in (1661) . 

Beyond the anticipated fact that the string leans back, the precise form of the 
relation we have obtained in 0651) is in fact directly associated with an important 
physical feature of motion in a strongly-coupled plasma. To see this, note first that 
the requirement that the Nambu-Goto action ([9]) describing a string embedded in the 
black hole metric ([6]) be real (i.e., that the string worldsheet be timelike) imposes a 
bound on classically realizable embeddings, which generalizes the familiar condition 
that a particle not be able to exceed the speed of light. The precise general bound 
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where 'd = Z{dX /dt,dX /dz). For a string that moves and stretches along a single 
direction x, as we are considering in the present paper, = or vr and fl67|) reduces 
to 

dxV , / , fdX 



The constraint fl68p . which is of course evident directly from fill I) , must be satisfied 
by the embedding function X{t,z) at all instants and all points along the string. In 
particular, at the physical string endpoint we must have 

< {1 + hms^) , (69) 

where hm = h{zm)- In terms of the external force [23], this reads 

h 



V < 



^/{h~^4JW+4J^) Vi< - AFV47rW[]4)(l + AFV47rW[]4) 

(70) 

where 

Vm = VK. = Vl - vr^T^^^ , (71) 

and [] denotes the expression within brackets in ([7]). 

In the unforced case we are presently considering, this implies the simple velocity 
bound 

V <Vm ■ (72) 

The appearance of (17T1) as a limiting velocity for quarks, mesons and other objects 
moving through a strongly-coupled thermal plasma is implicit in at least [l2l [211 122] , 
and was discussed more explicitly from various perspectives in [i3 | 1^ [36| [23l B5| H6] . 
The origin of the bound is easy to understand on the string theory side of the duality: 
a coordinate velocity v = v^. corresponds to a proper velocity V equal to that of light 
at the position of the string endpoint, z = Zm [12] • The interesting feature is that 
it is f , and not V , that corresponds to the gauge theory velocity. This yields then a 
nontrivial prediction of the AdS/CFT correspondence for strongly-coupled plasmas, 
which in turn might have interesting phenomenological consequences, such as the 
photon peak predicted in [SH] or the Cherenkov emission of mesons analyzed in [5U] . 
Now, for a heavy quark (consistent with fl43p ). the limiting velocity (17T]1 reads 

vm = l- In'T^i , (73) 
so we can finally recognize the physical significance of fl65|l : it states simply that 



V 

l-t; = Wm-f~l , (74) 

Vru. 
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which shows that the bound (172!) is equivalent to the standard relativistic constraint 
V <1. 

The preceding discussion motivates rewriting f l65p and f l66p in the form 

v = — (75) 

(a relation that is likely to be exact) and 




(76) 

which are more physically transparent and, to the order of our analysis, equivalent. 
Notice that the fact that 




VT^=\ , (77) 



implies that the predicted limiting velocity will be automatically incorporated into 
the kinematics of our finitely-massive quark. We will return to this point below. 
Notice also that (176|) and (165!) lead to 

~ 4t^4 4 



(1-^2)3/2 2(1-^2) 

which through (HUj) shows that our assumption (1501) is self-consistent, as promised. 

We comment in passing that (125]) and (1751) allow the embedding (1B2!) to be rewrit- 
ten in the form 



t{U,z) = U+ ^'Z^^ , (79) 



where all auxiliary variables have been eliminated and only the physical quark data 
remain. To avoid possible confusion, we should note that since the retarded time 
tr (where it is understood that x and v are to be evaluated) is defined just like in 
the unperturbed solution (l39l) . constant-t^ curves on the worldsheet are no longer null 
when finite-temperature contributions are taken into account. Given the way in which 
Mikhailov's vacuum solution ( II 7p (or its finite mass generalization (139 p ) incorporates 
information flow on the worldsheet at the speed of light, it would be natural to define 
a thermally-corrected retarded time by following the shifted null curves. This can 
be done, and amounts to a reparametrization of the worldsheet, but we will refrain 
from showing the result because it is not particularly illuminating at our level of 
approximation, even though it should presumably be crucial for the structure of the 
full finite-temperature solution. 
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We now finally proceed to extract the information on the evolution of the physical 
endpoint/quark. From fl28p we know that 



^2)3/2 



dtr = dtr ■ (80) 



Combining this with f l50|) . fl75|) and the fact that Vm is a constant, we see that, to 
lowest nontrivial order in our double expansion in small quantities, 

dv dv 
a = — = Vm-^ = Vma . 

Using (j751) . we deduce then that the quark/endpoint acceleration is given by 

1 W3 r-2 2 ^A^/^T^V^4^ .... 

« = -2^ ^ - = • (81) 

This is the finite-temperature (and unforced) analog of ( l40l) . and constitutes our main 
result. It shows that, as expected, the effect of the plasma is to slow down the quark, 
at an early-time rate that is naturally proportional to the h — 1 correction present in 
the Schwarzschild-AdS metric (|6]), and therefore differs markedly from that implied 
by the stationary/late-time result (Il])-(l2]). We will elaborate on this difference in the 
next section. 

The information contained in ( IHTl) can be alternatively cast in terms of a rate 
of momentum or energy loss for the quark, by imitating the procedure employed in 
Section [3?2l using f l75|) -f l80|) . the kinematic identity f l27|) . 

= fa , 

dtr 



can be rewritten in the form 




4^.4 3 2 

which is obviously equivalent to ( ISTl) . To interpret this in terms of momentum conser- 
vation, recall from Section [3l2] that, at zero temperature, the intrinsic quark momen- 
tum always takes the standard Lorentz- invariant form in terms of the auxiliary/tilde 
variables. Since these variables are defined at the AdS boundary {z = 0), where the 
presence of the black hole horizon is irrelevant, this same form must be preserved 
even in the finite-temperature context. In other words, the thermal corrections to the 
quark dispersion relation arise not from a modification to the standard form pq oc jv, 
but from the T-dependence present in the relations (I75l)- (176l) that connect the tilde 
and non-tilde variables. We thus identify 

mv 

Pg = m-fv = I ^ ^ (83) 



V 
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as the intrinsic quark momentum, and recognize (182|) as the statement (IMI) of mo- 
mentum conservation, with the external force F = and the rate of momentum loss 
given by 

dt 2{vl - v^) IQm^ivl -v^) ' ^ ' 

Actually, in the analysis of the previous paragraph we have ignored one subtlety. 
In Section IX^ we recalled from [23] the fact that the intrinsic quark momentum/energy 
arises from a worldsheet surface term. The same reference showed that, in the finite- 
temperature setting, this generally leads to an additional contribution from the upper 
string endpoint, i.e., the one at the spacetime horizon z = Zh- This is a thermal 
correction that was argued in [23] to encode the initial/boundary conditions of the 
quark -|- plasma system. For instance, in the case of a permanently static quark, 
the momentum surface term at the horizon vanishes, but the corresponding energy 
term gives the known —^/XT/2 correction to the quark rest mass. Precisely this same 
contribution appears in any quark configuration where the quark is assumed to have 
been at rest in the remote past, but the situation is different for other initial conditions 
(e.g., the case studied in [311] where the quark has constant velocity for all times [2B]). 
For concreteness, we will here assume the resting initial condition, which is in any case 
the most natural way to make contact between our setup and the phenomenological 
context, as we will discuss in the next section. Since —\/XT/2 is independent of time, 
at the level of the time-differentiated intrinsic quark momentum/energy appearing 
in the conservation laws ( IMl) and (l38l) . only the contribution of the lower [z = Zm) 
string endpoint matters, just as we implicitly assumed in the preceding paragraph. 
Notice that, with our assumed initial conditions, the undifferentiated expression fl83l) 
is correct as it stands. 

Bearing this point in mind, just like in Section [372] we can reinterpret the kinematic 
identity dj/dtr = ^^va as the energy conservation law ( !38|) . with -F = 0, 

Ei = '^^-7r—= r-, ^-^v^^ (85) 



the intrinsic quark energy, and 



dt 2{v^ — 16m2(f^ — 1>2) 



^6) 



as the rate of energy loss for the quark. Equations (!84l) and (l86l) are the results that 
we advertised in the Introduction. 

As noted below ( 177|) . we expect the quark kinematics (in the unforced case) to 
incorporate fl7ip as a limiting velocity, and this is indeed what is seen in fl83p and 
f l85|) : the intrinsic momentum/energy of the quark diverges as f — )■ v^M Notice also 
that, due to the connection with the auxiliary/tilde variables, for any type of initial 



^Of course, within our approximation scheme (|43p . we cannot get too close to this divergence. 
Our derivation does seem to suggest, however, that (ffS]). and therefore the divergence, is exact. 
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condition, the thermally-corrected dispersion relation for the quark always inherits 
the Lorentz structure 

iE,-Ei,y = {p,-p^f + m\ (87) 

where Eh and ph denote the corresponding horizon surface terms. This is in fact 
true even in the presence of external forcing, where, as studied in [23], the finite- 
temperature expressions for the quark intrinsic energy and momentum have to incor- 
porate the F-dependence seen already in the vacuum expressions and (15^ . This 
dependence is in particular responsible for a physical divergence at the critical value 
of the force associated with quark-antiquark pair creation [231 [191 [20] . Needless to say, 
it would be interesting to repeat the above derivation in the case with external force 
F 7^ 0, to deduce a thermally-corrected equation of motion that should presumably 
incorporate the full F-dependent bound ([70]) . Unfortunately, in the general case it is 
difficult to invert the expressions relating the tilde and untilde quantities. 

The equation of motion we have derived for the quark moving through the plasma, 
be it in the form ( IHT]) . ( [34]) or ( 138|) . admits an exact solution where v{t) is a cosine, 
of which only the linear portion of the evolution, 

^(t) = Vo-^^^vl-vit, (88) 

is within the regime consistent with our approximations. More precisely, to stay 
within our ultra- relativistic approximation fH5]) . we are limited to consider times 
such that the second term in fISS]) is much smaller than the first, i.e., 

2 16m^ , , 

t < tultra = , = , . (89) 



Notice also that, in view of f[79|l . to satisfy the small-disturbance condition (H3|l we 
must also have 

t < ttcmp = _ , , , (90) 



which is more restrictive than ( l89l) . because ( 143|) itself implies that ttcmp <S ^uitra- 



5 Phenomenological Estimates 

We will now try to make some quantitative inferences from our results in the previous 
section. For applications to the phenomenology of heavy ion collisions, we must choose 
values of the mass parameter based on the charm and bottom quark masses, 
m ^ 1.4,4.8 GeV. The issue of how best to translate between the SYM and QCD 
parameters has been discussed in [61] (see also [62]). Taking aqcD — 0.5 {qqcd — 
A/27r), Nc = 3 and Tqcd — 250 MeV, and employing the "obvious" prescription 
9ym = 9QCD {Qym^c ^ 19) and Tsym = Tqcd, we find from ([7]) that Zm/zh ~ 0.40 for 
charm and Zm/zh — 0.11 for bottom. If, on the other hand, one uses the "alternative" 
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scheme Qym^c — 5.5 (motivated in [61] through a rough matching of the AdS/CFT 
and lattice quark- antiquark potentials) and Tsym = ^^'^^'^Tqcd — 190 MeV (which 
follows from equating the energy densities of the two theories), then ([7j) leads to 
Zm/zh — 0.16 for charm and z^/zh — 0.046 for bottom. For definiteness, and knowing 
that the validity of our approximations would only improve if we used the bottom 
mass, we will take representative values in the neighborhood of the charm mass, 
Zm/zh = T^Tzm — 0.2-0.4, which is the same range studied numerically in |23]. We 
then have tt^T^z^ ~ 2 x 10~^ - 2 x 10~^, which is surely in compliance with fH^ . 

Next, we consider the velocities of the relevant quarks. For the mass values as- 
sumed above, the limiting velocity (17T]1 works out to Vm — 0.987-0.999. At RHIC, the 
presence of a parent hadron containing a single charm or bottom quark (i.e., a D or 
B meson, respectively) is inferred from the single electron resulting from its semilep- 
tonic decay Such electrons have been measured with transverse momenta up to 
Pt — 9-10 GeV, which, in the case of a charm decay electron, suggests that we may 
take 7 ^ 5 — 7 as representative values for the Lorentz factors. These correspond to 
quark velocities v ^ 0.980-0.990, which are indeed not too far from saturating the 
bound fl7^ . Using f l75]) . we thus have 1 — = 1 — [v/vmY ~ (2-4)xl0~^, in accord 
with our ultrarelativistic condition fHS]) . and (at least in the "alternative" setup) re- 
specting the ordering ( 149|) of the small parameters employed in the double expansion 
within which we have been able to derive analytic results. The numbers for LHC 
work out differently, since one expects heavy quarks there with transverse momenta 
up to ~ 100 GeV, but in any case, as noted in the Introduction, for extremely ener- 
getic quarks it will probably not be appropriate to carry out the entire energy loss 
calculation within a (weakly-curved) gravity dual. 

Having verified the pertinence of our approximation scheme, we are now set to 
explore the consequences of our results. The main question is how our early rate of 
momentum or energy loss f lH^ or f l86|) compares to the corresponding late/stationary 
rates ([2]) or ([T]), under conditions relevant to the quark gluon plasma as produced at 
RHIC. By taking the ratio between the two sets of expressions, we see that 



.2 



V dt lute _ \ dt ^latc _ 2!!iJL./T^r7p. r-^ - — ~2-6 r9n 

l^^J early l^^J early V ^ ^ ^"^ 

We see here that the initial damping proceeds at a rate that is a few times smaller 
than the asymptotic rate determined in [5l [6l [7]. This is consistent with the results 
obtained in [23] by studying the string evolution numerically, even though the analysis 
there was confined to nonrelativistic quark velocities. Given that the quark evolution 
at early times would be expected to be insensitive to the spatial extent of the plasma, 
it is interesting to note that a recent study [M] found that finite-size effects (relevant, 
by definition, at late times) also reduce the rate of energy loss as compared to the late- 
time infinite- medium results [3 |6l [7]. By way of contrast, the finite-mass correction 
to the stationary result of [3 [6] deduced in [ID] turned out to increase the rate of 
energy loss as compared to ([T]). 
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The appearance of in ( 186|) . (18^ and (IHT]) . which is sahently different from 
the T^-dependence seen in ([T]) and ([2]), and in our derivation clearly stems from the 
h — 1 deviation from pure AdS in the metric (jS]), is in consonance with an estimate 
for early-time energy loss presented in [2l], based on a heuristic argument involving 
saturation physics. By dimensional analysis in our conformally invariant theory, this 
scaling with temperature inevitably implies the dependence on zf^ seen in f l86|) . fl8^ 
and fl8ip . which in turn fixes not only the mass scaling but also the peculiar factor of 

Beyond the fact that the early-time damping rate obtained in this paper is smaller 
than its asymptotic/stationary counterpart, (!9T]) displays another interesting feature: 
the functional dependence on the quark velocity is substantially different in the two 
cases. Whereas in the late-time regime the friction coefficient defined through 



^ _ 1 dpg _ 1 dP.^d 

Pq dt Pq dt 



is momentum-independent [SlEj, 



/ilate = — , (93) 

2m 



in our early-time setup we find instead from flS^ and f l84l) that 

^"^''^ IQmH^^/vi^^ lQm% ' ^^^^ 

which given our restriction to the ultrarelativistic regime states that the friction co- 
efficient grows linearly with the quark momentum. This growth might seem counter- 
intuitive, but notice that flHT]) expresses the expected physical property that a faster 
quark is less damped by the plasmaQ When translated to the nuclear modification 
factor Raa{pt) (and similar observables) as in [151 HSl [13 HB], (IMj) would imply a 
PT-dependence that differs from that of both perturbative QCD and the late-time 
trailing string result ( p3|) . thus potentially serving as an experimental signature. 

In view of the above, a natural question is how far into the life of the quark-gluon 
plasma one could trust our early-time damping results (|5Ti) . (IMj) . and As 
discussed at the end of the previous section, the validity of our approximations (1431) 
and pHj) implies the temporal bounds fl89l) and fl90l) . of which the latter is the most 
restrictive. We thus know that we are limited to times much smaller than 

- . i^, (95) 

Given that I/ttT (for Tqcd — 250 MeV) corresponds to 0.25 fm/c under the "obvious" 
and 0.33 fm/c under the "alternative" prescription of [HI], we see that our limiting 



growing friction coefficient was also found in [69], for different physical reasons. 
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time (p5!) will be only of a few fm/c. This is an order of magnitude better than the 
time at which the numerical integration in |23] broke down (for the optimal scenario 
of a nonrelativistic quark), but is still, at best, of order the lifetime of the quark-gluon 
plasma, meaning that our results would only be justified in the early phase of the 
latter's evolution. 

For times of order ttcmp? we certainly expect the solution f l62|) to be modified 
by the appearance of and higher corrections neglected in the string equation of 
motion flH|) . As the string continues to evolve, we expect a dynamical worldsheet 
horizon (located at z = Zh in the remote past) to move down toward smaller values 
of z, and then (if the quark is undisturbed) move back up towards the spacetime 
horizon [23]. When the wavefront on the string crosses the worldsheet horizon, one 
would have to worry about the choice of boundary conditions there. Throughout this 
paper, we have employed the same condition as in [HI HSl [2D] , which describes waves 
in the gluonic field that are required to propagate outward from the quark. This 
captures the causal setup that is most natural for a quark in vacuum, and is therefore 
also appropriate as long as we consider only small perturbations about f lT7|) . as we 
have done here given that our goal from the outset was to find the first correction 
due to the thermal medium. When the evolution progresses further, however, it is 
to be expected that, as the medium is disturbed by the quark, it radiates gluonic 
waves back towards the source. Needless to say, determining and implementing the 
precise boundary condition in a fully dynamical setup is bound to be a challenging 
problem, which would seem to require in particular knowledge of the exact solution of 
the nonlinear thermal equation (H2l) . The net result of the additional thermal effects 
described in the present paragraph should be to yield a damping coefficient /i(t) that 
smoothly interpolates between the early-time result flMl) and the late-time form (p3|) . 

Beyond the limitation to early times, it should be borne in mind that the setup in 
this paper, as well as in O |6l [7] , involves an infinite static thermal medium, whereas 
the real- world quark-gluon plasma is of finite extent, and is exploding and cooling 
rapidly. Starting with [51], the literature contains various examples of gravity duals 
for expanding plasmas (for reviews, see, e.g., [66]), but given the reliance of our ap- 
proach on the exact zero-temperature embedding (ITTj) derived in [H] by exploiting 
the symmetries of the pure AdS background, at present there seems to be unfortu- 
nately little hope of being able to extract from them analytic results for the early-time 
thermal damping rate (see, however, [ST] for related results). 

Another issue is the choice of initial conditions for the quark. We wish to model 
a situation where the quark is created within the plasma and is only starting to feel 
its effects. To account for this, on the gravity side we have been considering a string 
that is disturbed only in a region close to its endpoint, with the disturbance thereafter 
traveling deeper into the bulk, but not having had time yet to get close to the black 
hole horizon. The cleanest scenario involves a string that is originally at rest, and 
whose endpoint is then accelerated via an external force for some short period of 
time, after which it is released. This is precisely the situation that was examined 
numerically in [23]. The period of external forcing can be viewed as modeling the 
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formation time for our energetic quark within the plasma, and it merely serves to set 
some initial conditions. Both from general physical principles and from the numerical 
results of |23], we know that, if this period is short enough, the quark will behave as 
in vacuum, and so the corresponding evolution of the string will be essentially as in 
pure AdS. The question we have been addressing in this paper is what happens after 
the quark is released and slowly begins to feel the effect of the thermal medium. 

Another possibility that might come to mind would be to imagine a quark traveling 
at constant velocity in vacuum, which then suddenly encounters the thermal medium. 
A similar setup was explored in [211 127]. In this case, in the gravity description we 
would start with a vertical string translating uniformly in pure AdS and then in some 
fashion 'switch on' the h — 1 metric deformation leading us to the black hole geometry. 
However, in such a situation, the portion of the string that first feels the thermal 
effects would clearly be the one closest to the spacetime horizon ai z = Zh, where our 
condition f H3|) is violated. It is easy to see that the vertical string traveling at constant 
velocity is in fact a solution of the full thermal equation of motion f H2|) . so that is 
not an issue. The problem, known already from O |6], is that this string worldsheet 
becomes spacelike above the critical radial position = {1 — v'^Y^^Zh, and is therefore 
unphysical. So, in this scenario where the quark is suddenly introduced into the 
plasma, there will necessarily be a wave running downward along the string to ensure 
that a physical boundary condition is achieved at the emerging worldsheet horizon. 
This is the same type of effect that we argued above to be important to eventually 
lead to the late-time/stationary rate ([T]), and for the same reason, it is beyond the 
reach of our approach. In MSYM language, the initial quark in this scenario has 
its full vacuum dressing (i.e., the gluonic cloud surrounding it is undeformed out to 
arbitrarily large distances), and so does not closely mimic the experimental situation. 

Returning to the scenario where the string is initially static and is then forced over 
a finite period of time, we should emphasize that our use of the vacuum configuration 
as a starting point does not restrict our analysis to situations where the string remains 
nearly vertical (which would in turn require the forcing not to be too violent), because 
Mikhailov's solution (fTTj) is valid for arbitrary trajectories of the quark. The vertical 
string profile is recovered in the case where the quark moves at constant velocity, 
but an exact description is equally available in situations where the quark is forced 
(even violently), such as uniform proper acceleration (where the solution was found 
independently in [23 EH] but was later shown [S3] to be equivalent to ( ITT]) ) or uniform 
circular motion (where, again, the solution derived by a different route in [69] actually 
coincides with Mikhailov's). 

Now, in the preceding discussion we have for simplicity been considering the evolu- 
tion of an isolated quark as it travels through the plasma, but since the quark should 
be initially created with its corresponding antiquark, we could expect the latter to 
have some influence on the former. This concern has been addressed in [23], follow- 
ing [5] (related, more recent, work may be found in [70]). The gravity description 
involves two string endpoints on the D7-branes, which separate from one another. As 
explained in [23], there are in fact two cases to consider. A quark-antiquark pair in 
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the adjoint representation, such as would be created by a gluon, corresponds to a A- 
shaped string with its apex located at z = Zh, that is initially of zero width and then 
proceeds to open up as the string endpoints move away from each other. In this case, 
since the two halves of the string — corresponding to the quark and antiquark field 
configurations — only meet at the black hole horizon, they do not interact in finite 
time (at large- iVc), and we are back to the isolated quark scenario described above 
(albeit with some differences introduced by the details of the forcing stage). This 
is then the most realistic interpretation of the initial conditions naturally associated 
with our approach. 

The other case refers to a quark-antiquark pair in the singlet representation, such 
as would be created by a photon, and corresponds to a fl-shaped string that is initially 
of zero length, with its midpoint moving deeper into the bulk as the two endpoints 
separate. In this case, it was found numerically in [23] that the evolution is essen- 
tially as in vacuum up to the point where the quark and antiquark lose contact by 
moving beyond their (velocity-dependent) screening length [171 [SH HH], after which 
their evolution is already controlled by the late-time damping ([1]). Our results (or 
their generalization to finite external forcing) would then be expected to be relevant 
only for the initial stage, prior to screening. In this connection we should note that, 
even though the general vacuum solution (fTTj) is nominally associated with an embed- 
ding that has a single endpoint on the D7-branes, and therefore describes an isolated 
quark, it is in fact also capable of describing at least certain types of fl-shaped profiles, 
with two endpoints on the D7-branes, such as the example already mentioned two 
paragraphs above, where the quark and antiquark separate back-to-back with con- 
stant proper acceleration. Again, the period of external forcing we have envisioned 
in our above discussion on initial conditions can be regarded as modeling the stage 
of formation of the quark, and, in particular, the influence of the antiquark. 

Altogether, then, we believe the results obtained in this paper shed some addi- 
tional light on the nature of early-time energy loss in a scenario that is not too distant 
from the one that is of phenomenological interest, but of course, additional work is 
needed to achieve further realism. 
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